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Summary 

A technique called relative control effectiveness 
is developed to give a quantitative measure of the 
effectiveness of each of several controls upon the 
fundamental modes of a linear dynamic system. By 
postulating linear combinations of controls, a set 
of pseudo controls is formed. The relative control 
effectiveness measure is then used to formulate a set 
of optimization problems which yield the proportions 
of the controls to be assigned to each of the pseudo 
controls with the objective that each pseudo control 
strongly affects selected fundamental modes of the 
system while only weakly affecting the remaining 
modes. 

The procedure is applied to the linearized lateral 
dynamics of a high-performance fighter airplane to 
obtain (1) a pseudo control which primarily affects 
the roll and spiral modes and (2) a second pseudo 
control which primarily affects the Dutch roll mode. 
Transient responses of the linear airplane model to 
lateral and directional aerodynamic controls and to 
the pseudo controls are presented^ An example of the 
application of the procedure to airplane control with 
thrust vectoring is included. 

Introduction 

The usual definition of controllability of a linear 
system may be stated as follows: 

A linear system is controllable by a given 
control variable if that variable can excite each 
fundamental mode of the system. (See refs. 1 
and 2.) 

If this criterion is true for a given linear system, it is 
mathematically possible to devise a control scheme 
which influences all of its fundamental modes. How- 
ever, this criterion does not indicate the degree of 
difficulty which would be experienced in generating 
such a control scheme. For example, a mode may 
be only slightly affected by reasonable, but limited, 
control deflections. Such a system would be “con- 
trollable” by the above definition but would be un- 
controllable in practice. 

What is needed is a quantitative measure of con- 
trollability which would indicate the “degree of con- 
trollability” of the system to control deflections. This 
measure could be used to compare the control effec- 
tiveness of multiple control variables upon each of 
the modes of the system and to suggest a preferred 
control structure; i.e., which control variables should 
be used to control each mode. Multiple control vari- 
ables could be mixed in linear combinations, with the 
proportions of the mixings being determined by plac- 
ing emphasis upon the quantitative controllability of 


specific modes. This approach could be useful in de- 
termining control-mixing schedules in the presence of 
changing operating conditions. The technique could 
also be used to reduce the number of control channels 
in the case of multiple redundant controls (for exam- 
ple, on airplane control using multiple aerodynamic 
controls and thrust vectoring). 

A relative control effectiveness measure was de- 
veloped and is described in this report. A method 
of applying this measure to obtain control-mixing 
combinations to provide control channels which can 
(more or less) independently affect the fundamental 
modes of a linear system is presented. The method 
is applied to an example of the lateral control of 
a modern high-performance fighter airplane in low- 
speed flight. This report also includes an example 
for thrust- vectoring nozzles. 

Symbols 

A state variable coefficient matrix (n x n) 

a c i real part of ith complex eigenvalue of A 

a T i ith real eigenvalue of A 

B control variable coefficient matrix (n X m) 

b c i imaginary part of ith complex eigenvalue of 
A 

C m x k matrix of control-mixing vectors (c j) 

C e fi relative control effectiveness matrix (num- 
ber of modes x m) 

Cj jth control-mixing vector of dimension m 

J quadratic cost function for maximization 
problem 

k number of pseudo control variables 
M similarity transformation matrix ( n x n) 
m number of control variables 

n number of state variables 

p body axis roll rate, deg/sec 

r body axis yaw rate, deg/sec 

t time, seconds 

TAS true airspeed, knots 

u vector of control variables (iq) of dimension 

TO 

u vector of scaled control variables (u,) of 

dimension to 

v vector of pseudo control variables ( Vj ) of 
dimension k 



w mxm matrix of control variable authorities, 
W = diag(iOi) 

W symmetric matrix for maximization prob- 
lem (m x m) 

Wi maximum value of control variable u t 

x vector of state variables (27) of dimension n 

y vector of transformed state variables in 
modal coordinates (y t ) of dimension n 

a angle of attack, degrees 

(3 sideslip angle, degrees 

r transformed system control variable coeffi- 

cient matrix (n x m) 

7, column vector formed of elements of the ith 
row of matrix T 

7 ij element of matrix T in the ith row and jth 

column 

6 a aileron deflection, degrees (positive, right 

aileron up) 

6 d differential horizontal-tail deflection, de- 

grees (positive, right tail trailing edge up) 

<5 r rudder deflection, degrees (positive, trailing 

edge left) 

<5u,roii rolling thrust vector, degrees (positive, right 
nozzle up) 

(5., yaw yawing thrust vector, degrees (positive, 
nozzles left) 

6, iP phase angles 

A transformed system state variable coeffi- 
cient matrix in block diagonal form (n x n ) 

A Lagrange multiplier 

/ i c i eigenvector of dimension n corresponding to 

ith complex eigenvalue of A 

Hri eigenvector of dimension n corresponding to 
ith real eigenvalue of A 

52 summation over desired modes 

D 

52 summation over nondesired modes 

ND 

r integration variable 

4> bank angle, degrees (positive, right wing 

down) 

A dot over a variable (e.g., x) denotes the deriva- 
tive with respect to time. 


Control Effectiveness Measure 

In the analysis and design of multivariable con- 
trol systems, it is customary to express the small- 
perturbation system dynamics in the form of a set of 
coupled linear first-order equations in the following 
matrix representation: 

x = Ax + Bu (1) 


where 

x vector of state variables of dimension n 
u vector of control variables of dimension m 
A n x n matrix of state variable coefficients 
B n x m matrix of control variable coefficients 

The state variable coefficient matrix, A, deter- 
mines the characteristics of the fundamental modes 
of the system. The control variable coefficient ma- 
trix, B, specifies the influence of the control vari- 
ables on the state variables. The matrix B does 
not directly indicate the controllability of the system. 
Moreover, a column of matrix B with large elements 
does not necessarily indicate that the system is in- 
fluenced by the corresponding control variable any 
more than by any other control variable. 

The controllability criteria of references 1 and 2 
can be applied to equation (1). These criteria indi- 
cate whether it is mathematically possible to influ- 
ence all the modes of the system by manipulation of 
the control variables. If not, the criteria identify the 
modes which cannot be altered by any means. Un- 
fortunately, an ill-conditioned system may satisfy the 
controllability criteria but possess one or more modes 
which are weakly affected by all the control variables. 
Attempts to excite or modify such modes without un- 
duly upsetting the remainder of the system may be 
futile in practice. This report presents an alternative 
method of measuring controllability which includes 
quantitative information in addition to the qualita- 
tive information of the conventional controllability 
criteria. 

An indication of the influence of the control vari- 
ables upon the system modes may be obtained by 
transforming the system equations to a block diago- 
nal form 

y = Ay + Tu (2) 


A = M -1 AM : 


“rl 


a r2 


a c\ b c 1 
~ 1 a c 1 


( 3 ) 
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where 


r = M BW 

(4) 

x = My 

(5) 

u = Wu 

(6) 

W = diag(iUi) 

(7) 


The similarity transformation matrix, M, is cho- 
sen such that it is analogous to the modal matrix 
used in references 1 and 2, except that it is a real ma- 
trix. The matrix M, which is sometimes called the 
right eigenvalues of A, has real coefficients because of 
the different method of handling the complex modes. 
This choice results in the transformed state variables 
y (the modal coordinates) and the coefficient matri- 
ces A and T also being real. 

The similarity transformation matrix, M, is as- 
sembled from the eigenvectors of A. For each real 
eigenvalue of A, one column of M is set equal to the 
corresponding eigenvector as is done in constructing 
the modal matrix. For each complex conjugate pair 
of eigenvalues of A, one column of M is set equal to 
the real part of the corresponding eigenvectors, and 
the adjacent column is set equal to the imaginary 
part of either eigenvector. In mathematical notation, 
this is 


M=[/i rX , ..., Re(#i cl ), Im(/i cl ), . . .] (8) 

where fi ri is the eigenvector corresponding to the 
zth real eigenvalue and fi c , is a complex eigenvector 
corresponding to the zth complex conjugate pair of 
eigenvalues. 

The matrix W has the maximum value (author- 
ity) of each of the control variables (wi = max(zq)) 
as its diagonal elements. The transformed system is 
driven by the scaled control variables v,i which have 
unity ranges (±1). 

In the transformed state coefficient matrix of 
equation (3), a ri is the zth real eigenvalue, and a cl 
and b c i are the real and imaginary parts, respectively, 
of the zth complex pair of eigenvalues. The elements 
of y, which are commonly referred to as modal 
coordinates, are indications of the excitation of the 
modes. The influences of each control variable upon 
the fundamental modes are given by the elements 
of the transformed control coefficient matrix, T, of 
equation (4). The absolute sizes of the elements of 
r have little significance in themselves. However, 
they are useful in assessing the relative influence 
of the various control variables upon each of the 
modes. Since scaled control variables are being used, 


the elements of T give an indication of the control 
effectivenesses upon the fundamental modes, which 
is consistent with respect to the authorities of the 
control variables. 

For a real mode, one may compare the elements 
of the corresponding row of the matrix T to get an 
immediate comparison of the relative effects of the 
control variables upon that mode. For a complex 
mode, there is a pair of transformed state variables, 
?/i and y l+ x , which are coupled. Therefore, there are 
two rows of the matrix T to be considered. Consider 
the following second-order example: 



where the subscripts on a and b have been omitted 
for brevity. The general solution of this system may 
be written as 



where 


+ 2/2(°) 

at j cos(6t - 6) ) 
e ( sin(6t — 0) J 

■ 

+ Ev^ 

j 

+ K;) 2 



cos[fc(t — r) — W] ] 
sin[6(t — t) — 'J] J 

| “j( r ) dr 



(10) 

6 = tan -1 

/ 1/2(0) \ 
U(0)J 

(ii) 

$ = tan - 

V7i j7 

(12) 


It is evident that the control variables influence 
the two transformed state variables in equivalent 
fashions and that this influence is proportional to 

yjllhj) 2 + il2j) 2 ■ 

A relative control effectiveness matrix which in- 
corporates these concepts can now be defined. For 
each real mode, the relative control effectiveness of 
the jth control variable upon the zth first-order mode 
is taken to be the magnitude of the (*, j) element of 
r divided by the Euclidean norm of the zth row of T, 


C'eff (b J ) 




(13) 


For each complex mode, the relative control effective- 
ness of the jth control upon the (z,z + l)th second- 
order mode is taken to be the square root of the sum 
of the squares of the (z, j) and (z + 1, j) elements of T 
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divided by the Euclidean norm of the fth and (i-f-l)th 
rows of r combined, 




l + ('7i+lj) 2 
E f(7.o) 2 + ('7»+i,i) 2 | 

3 


(14) 


The elements of the control coefficient matrix, T, 
provide a somewhat direct indication of the effective- 
ness of each of the control variables in influencing 
each of the fundamental modes of the system. How- 
ever, the magnitudes of the elements are dependent 
upon the lengths of those eigenvectors used to trans- 
form the system into block diagonal form. The el- 
ements of the relative control effectiveness matrix, 
C e ff, are limited to positive values between zero and 
unity and are independent of the eigenvector lengths. 
However, they do not provide any absolute informa- 
tion about the effectiveness of the controls. 


Specification of Control Crossfeeds 

Systems having multiple control effectors provide 
numerous paths by which control law schemes may 
direct signals from the operator controls and the sys- 
tem output variables to the control variables. Mul- 
tivariable design techniques, such as those imple- 
mented in reference 3, generate a large number of 
control gains. Only a few of these may be significant. 
Attempts to systematically eliminate the ineffectual 
gains have been made in references 4 through 6. 

The approach taken in the present study is first 
to find control-mixing combinations (control cross- 
feeds). A number of “pseudo” controls are formed. 
The activation of each of the pseudo controls com- 
mands proportional deflections of the system con- 
trols. Once a suitable set of pseudo controls is estab- 
lished, the design of control laws can be conducted 
by means of conventional techniques. 

For systems which have multiple, redundant con- 
trols, this procedure reduces the number of controls 
which must be dealt with in designing control laws. 
The control system structure is simplified, and the 
number of control gains is reduced. If the pseudo 
controls are formulated so that each has its effect 
concentrated upon a single selected mode, then the 
control law design problem decouples into a series of 
first- and second-order problems. In general, the de- 
coupling will not be complete but may be adequate to 
allow the utilization of a simplified control structure. 
A technique for calculating the required relationships 
between the pseudo controls and the system control 
variables is developed next. 

Define a set of pseudo control variables, Vj, which 
are linearly related to the scaled control variables, ft,, 


through the matrix C such that 
u = Cv = c j 


(15) 


where the columns of C have unity length; i.e. , 
c Jcj = 1. Combining equations (15) and (2) yields 
the following form for the transformed system: 


A y + E 




T 

ll 


y In 


(16) 


where qr J is the fth row of I\ The vectors c j are 
control-mixing combination vectors which specify the 
proportions in which the system control effectors are 
actuated. 

For each pseudo control, v v a control-mixing 
vector, Cj, is found. The effect of the pseudo control 
upon one or more modes is maximal while the effect 
on other modes is minimal. To accomplish this, the 
following maximization problem is solved. Find the 
control-mixing vector, c which maximizes 

J = cj ^ c > ~ c 3 E ll1 ' Cj ( 17 ) 

D ND 

subject to the constraint 

cjcj = l (18) 

where the first term of J includes the mode(s) to 
be controlled by Vj (desired), and the second term 
of J includes the modes to be unaffected by Vj (not 
desired). 

This constrained maximization problem is solved 
by introducing a Lagrange multiplier, A, and finding 
the combination of c j and A which maximizes the 
following expression: 



A(cJ c j - !) (19) 


Taking the derivatives of this expression with re- 
spect to Cj and A yields necessary conditions for the 
solution 


i^=Wc i -Ac J -=0 (20) 

2 dcj 3 
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where 


dJ T 



< p i 


- - 2.13 

0 

2.19 

- 12.7 ' 


< p 

ax =c ’ c ‘~ 1=0 

(21) 

d 

— 

1 <t> \ 

> = 

1 

0 

0.176 

0 

< 

I <t> 



dt 

r 1 


0.0646 

0 

- 0.559 

1.44 


r 




UJ 


. 0.174 

0.114 

- 0.985 

- 0 . 160 . 


. 0 . 


w = £ 7,7 t T - ^ IrlJ (22) 

D ND 

Equations (20) and (21) are in the form of an 
eigenvector decomposition problem. The control 
mixing vector, c j, is a normalized eigenvector of ma- 
trix W, with A being its associated eigenvalue. Since 
W is real and symmetric (eq. (22)), its eigenvalues 
and eigenvectors are real and orthogonal (ref. 2). By 
combining equations (19) through (22), the value of 
J is found to be equal to the eigenvalue, as follows: 

J = cJ\Vcj = cj Acj = A (23) 



' 4.38 

1.10 

4.09 ' 

1 . . . 


0 

0 

0 

i 

r a 

+ 

- 0.212 

- 1.20 

0.222 




.- 0.00155 

0.0330 

- 0 . 0103 . 

t 


( 24 ) 


where 

p roll rate, deg/sec 

<f> bank angle, degrees 

r yaw rate, deg/sec 

j3 sideslip angle, degrees 


Therefore, the maximizing solution for Cj is the one 
for which the corresponding eigenvalue, A, has the 
maximum positive value. 

The technique is applied by partitioning the sys- 
tem modes into groups. A pseudo control which has 
its principal effect upon each group is to be found. 
The control-mixing matrix for each pseudo control 
is then found by solving the above eigenvector prob- 
lem. After this is accomplished, any of the usual 
multivariable control design techniques may be used 
to complete the design of the control system. 

Throughout the above analysis, it is implicitly 
assumed that the fundamental modes were equally 
important. If this is not the case, the procedure could 
be extended in a straightforward way to apply a set 
of weights to differentiate the fundamental modes. 


Example Applications 


6 a aileron deflection, degrees 

6 r rudder deflection, degrees 

S D differential horizontal-tail deflection, 

degrees 

The block diagonal (modal) form of this system 
(cf. eqs. (2) through (8)) is 


■- 

-1.37 

0 

0 


0 ' 


0 

-0.101 

0 


0 


0 

0 

-0.693 1.56 

> 

0 

0 

-1.56 

-0.693 J 


’ 156 

-124 

no - 


' 6J 40 ) 
<5 r /30 \ 

{ W20 J 

+ 

43.0 

38.7 

-120 

-3.67 

51.4 

21.4 

< 


. 147 

212 

23.3. 



y 


(25) 


The concept of relative control effectiveness and 
a technique to select control-mixing combinations for 
generating those pseudo control variables which prin- 
cipally affect selected modes of a system have been 
developed. This section applies these developments 
to two example flight control problems. 


where the similarity transform matrix is 


' 0.805 -0.116 -0.654 

-0.591 0.988 0.437 

0.017 0.092 0.013 

.-0.046 0.034 0.144 


0.548' 
0.233 
-0.107 
0 . 


(26) 


Lateral-Directional Airplane Controls 

Consider the lateral control of a modern high- 
performance fighter airplane in level flight at a true 
airspeed of 165 knots with an angle of attack of 
10°. This flight condition approximates the normal 
landing approach speed for the example airplane. 
The lateral perturbation model for this example (cf. 
eq. (1)) is 


The authorities of the aileron, rudder, and differen- 
tial tail are 40°, 30°, and 20°, respectively. The first 
real mode of this system is the roll mode, the second 
real mode is the spiral mode, and the complex mode 
is the Dutch roll mode. 

The relative control effectiveness matrix, C e fj, is 
calculated by applying equation (13) to the first two 
rows of the control coefficient matrix, T, appearing 
in the second term of equation (25) and by applying 
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equation 


C e ff = 


(14) to the second two rows as follows: 



Sr 

S D 


0.684 

0.544 

0.486 

Roll 

0.312 

0.874 

0.373 

Spiral 

0.577 

0.808 

0.120 

Dutch roll 


(27) 


Examination of this matrix indicates that the effects 
of the three controls upon the roll mode are com- 
parable. The rudder is dominant in its influence on 
the spiral and Dutch roll modes, and the ailerons 
strongly affect the Dutch roll mode. Unless the con- 
trols are carefully proportioned, those control com- 
mands generated by either cockpit controller move- 
ments or by feedback control mechanisms will excite 
all the modes of the system and result in potentially 
undesirable intermode couplings. 

Two pseudo control variables were calculated for 
this system. It was intended that the first would have 
its principal effect upon the roll and spiral modes, 
and the second, upon the Dutch roll mode. The 
matrices for the eigenvalue decomposition problems, 
W, are calculated according to equation (22). For the 
roll and spiral pseudo control variable, v h the first 
term of equation (22) includes the first two rows of 
r, and the second term includes the second two rows. 
The calculated W matrix follows: 


W = 


3070 
-55 400 
15 200 


-55 400 
-15 300 
-24 700 


15 200 
-24 700 
13 900 


(28) 


The eigenvalues of equation (28) are 65 100, —0.1, 
and —63 400. The control-mixing vector for the 
roll and spiral modes is the eigenvector of equa- 
tion (28) corresponding to the largest positive eigen- 
value (namely, 65 100). 

For the Dutch roll pseudo control variable, r> 2 , the 
terms of equation (22) are interchanged for this ex- 
ample. The matrix W is equal to the negative of 
equation (28). Therefore, the control-mixing vector 
is the eigenvector of equation (28) corresponding to 
the most negative eigenvalue (namely, —63 400). It 
is fortuitous that both control-mixing vectors can be 
found by solving a single eigenvalue decomposition 
problem for this example. In general, a separate de- 
composition problem must be solved for each pseudo 
control variable. 

The eigenvectors found above are arranged to 
form the control-mixing matrix (cf. eq. (15)) as 
follows: 



25.9 

24.7' 

\t>r = 

-17.8 

23.3 

l I 

9.56 

2.54 



(29) 


Note that the control scale factors have been in- 
cluded in equation (29), so that 6 a , 6 r , and S D are 
in terms of degrees of deflection, whereas the pseudo 
control variables are dimensionless with ranges of 
approximately ±1. The numerical procedure used 
to calculate the eigenvectors of equation (28) might 
cause sign inversions in any of the columns of equa- 
tion (29). It may be necessary to change the signs of 
the columns in order to obtain the desired phasing of 
the controls. 

Use of equation (25) and the control-mixing ma- 
trix of equation (29) gives the relative control effec- 
tiveness matrix for the airplane, 



v 2 


0.998 

0.062' 

Roll 

0.899 

0.438 

Spiral 

0.161 

0.987 

Dutch roll 


(30) 


The roll and spiral modes are now principally affected 
by Vi with small coupling of the Dutch roll mode. 
This pseudo control causes the aerodynamic controls 
to be activated in a coordinated fashion. A positive 
value of v x yields (1) aileron and differential tail, left 
surfaces trailing edge down and (2) rudder, trailing 
edge right. The Dutch roll mode is principally af- 
fected by i> 2 with some coupling of the spiral mode. 
This pseudo control causes crossed controls. A posi- 
tive value of v 2 yields (1) ailerons and differential tail, 
left surfaces trailing edge down and (2) rudder, trail- 
ing edge left. The partial decoupling of the system 
modes improves the conditioning of the system for 
manual control and for the design of feedback con- 
trol laws. It has also reduced the number of control 
variables which must be considered in controlling the 
airplane. 

Transient responses were obtained for control step 
inputs for the two pseudo controls and for lateral 
(ailerons with differential horizontal tail) and direc- 
tional (rudder) aerodynamic controls. The control 
effector deflections used are given in table I. These 
deflections are 10 percent of the available values. The 
time histories of the airplane state variables are pre- 
sented in figures 1 and 2. These data were obtained 
from analog computer simulations of the linearized 
dynamics of equation (24). In each case, the controls 
were deflected for approximately 6 seconds and then 
returned to their neutral positions. 

Figure 1(a) shows the airplane open-loop response 
to simultaneous deflections of the ailerons and differ- 
ential tail with a 2-to-l gearing ratio. The presence 
of the Dutch roll mode appears as a damped oscil- 
lation in the roll rate and sideslip responses. The 
slow exponential character of the spiral mode can be 
seen in all the state variables. Figure 1(b) shows the 
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airplane response to rudder deflection. The airplane 
reaches a fairly steady roll rate and sideslip after an 
initial overshoot. The spiral mode is also present in 
all the state variables. 

Figure 2(a) shows the airplane response to the 
pseudo control variable rq, which was selected for 
maximum influence on the roll and spiral modes. The 
responses are quite different from those of figure 1(a). 
The presence of the Dutch roll mode is not apparent, 
the roll rate achieves a significantly greater magni- 
tude, the sideslip smoothly follows the bank angle, 
and the yaw rate responds immediately in the de- 
sired direction. 

Figure 2(b) shows the airplane response to the 
pseudo control v 2 , which was selected for maximum 
influence on the Dutch roll mode. These responses 
are different from those of figure 1(b). The roll rate 
is seen to settle to a smaller value, and the sideslip 
angle is increased and achieves a steady value. 

Aerodynamic and Thrust- Vectoring Controls 

In this example, the mathematical model of 
the airplane was modified to include ideal thrust- 
vectoring nozzles at the extreme rear end of the 
airplane to provide additional rolling- and yawing- 
moment control capability. Each nozzle was assumed 
to be capable of changing the direction of thrust 15° 
in any direction. Rolling moments are produced by 
vectoring one nozzle upward and the other nozzle 
downward for an available deflection (5 V)ro n) of 30°. 
Yawing moments are produced by vectoring the noz- 
zles simultaneously to the side for an available de- 
flection (<5 v , yaw ) of 15°. 

The airplane was trimmed for steady level flight at 
an angle of attack of 15° and an airspeed of 133 knots. 
This airspeed is approximately 32 knots less than the 
normal landing approach speed for flaps up, and the 
angle of attack corresponds to the maximum pitch 
attitude allowed at touchdown to avoid tail skid. 

The control effectiveness technique was applied to 
this model to find control-mixing combinations which 
maximize the responses of the roll and spiral modes 
with minimum excitation of the Dutch roll mode. 
The technique was applied to the case in which both 
the aerodynamic and the thrust-vectoring controls 
are active and to the case in which only the aero- 
dynamic controls are active. The calculated control- 
mixing matrices and relative control effectiveness ma- 
trices for this example are given in tables II and III. 

The value of the pseudo control variable which 
primarily affects the roll and spiral modes, was se- 
lected so that the deflections of all the airplane con- 
trols would be less than or equal to 15 percent of 
their authorities. These deflections are presented in 


table IV. Also included for reference is a baseline 
example which has the ailerons and differential tail 
deflected to 15 percent of their authorities. Tran- 
sient responses of the airplane (linear model) were 
obtained for step inputs of the lateral controls with 
values as given in the table. These step responses are 
presented in figure 3. 

The baseline example (fig. 3(a)) is objectionable 
for several reasons: (1) the oscillatory nature of the 
roll rate, which causes the bank angle to “ratchet,” 

(2) the initial negative response of the yaw rate, 

(3) the rapid rise of the sideslip angle, and (4) the 
small constant value reached by the bank angle, 
which implies the necessity of holding a substantial 
lateral-control input in order to maintain a desired 
bank angle. This behavior of the baseline example 
can be attributed to the adverse yaw characteristic 
of the ailerons which excites the Dutch roll mode. 

Smooth increases in bank angle with well-behaved 
yaw rate and reduced sideslip are obtained by using 
mixed aerodynamic controls (fig. 3(b)) and mixed 
aerodynamic and thrust- vectoring controls (fig. 3(c)). 
The airplane turn coordination has been improved. 
The limiting control for these cases is the rudder. 
The case of figure 3(c) uses a relatively large amount 
of the yawing thrust vector to supplement the rudder. 
The rolling thrust vector is ineffective because of the 
small distance between the engine nozzles. Since it is 
ineffective, the deflection of the rolling thrust vector 
is made small in comparison with the more effective 
controls. 

The addition of thrust-vectoring controls (to the 
airplane of this example) has increased the roll rate 
capabilities of the airplane by approximately 60 per- 
cent with only a small increase in sideslip while us- 
ing moderate deflection angles. The improvement 
of the airplane response was achieved by appropri- 
ately gearing the controls without using any feedback 
loops. 

Concluding Remarks 

A method for selecting mixing combinations 
(crossfeeds) for the control variables of a linear sys- 
tem has been developed. The mixing combinations 
are chosen so that the resulting control channels have 
their principal influences on selected fundamental 
modes of the system. In the ideal case, each fun- 
damental mode would be controlled by one control 
channel independent of the others. The control de- 
sign problem would then reduce to a series of first- 
and second-order design problems. This ideal decou- 
pling is not possible in general. A series of algebraic 
maximization problems is used to find the control- 
mixing combinations which maximize the effects of 
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the control channels on selected modes while simul- 
taneously minimizing the effects on the remaining 
modes. 

The technique allows a reduction in the number 
of channels in the case of multiple, redundant con- 
trols. The most effective control effectors are selected 
in the proper combination for each operating con- 
dition. Given a judicious selection of modes to be 
associated with a number of pseudo controls, the re- 
sultant control-mixing combinations may transform 
the system into one which is better conditioned for 
the subsequent application of feedback control law 
design procedures. 

A quantitative measure of relative controllability 
is a product of the procedure. This measure can be 
used to assess the influence of each control variable 
on the modes of the system as well as on the relative 
effectiveness of the control-mixing combinations. 

The method was applied to the lateral and di- 
rectional control of an example airplane trimmed 
in level flight at an angle of attack of 10°. Two 
pseudo control variables were calculated for this sys- 
tem. The first has its principal effect on the roll 
and spiral modes, and the second, on the Dutch roll 
mode. Use of these pseudo control variables was 
found to eliminate the oscillations present in the roll 
rate for a step lateral-control input and to improve 
the sideslip response with reduced rolling motions for 
a step directional-control input. 


The airplane model was modified to include 
thrust-vectoring nozzles which provide rolling and 
yawing control moments for use in very slow flight 
at moderate angles of attack. Yawing moments gen- 
erated by thrust vectoring significantly improve the 
roll rate capability of the airplane. The rolling mo- 
ment generated by thrust vectoring was found to be 
ineffective for the example studied. 

NASA Langley Research Center 
Hampton, VA 23665 
December 11, 1984 
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TABLE I. CONTROL DEFLECTIONS: LATERAL-DIRECTIONAL 
[a = 10°; TAS = 165 knots] 


Control 

variable 

(authority) 

Control deflection, degrees 

Baseline system 

Pseudo controls 

Lateral 
(fig. 1(a)) 

Directional 

( fi g- 1(b)) 

(fig. 2(a)) 

V2 

(fig- 2(b)) 

Sa (40°) 

*4 


*4 

3.18 

6r (30°) 


*3 

-2.75 

*3 

6 d (20°) 

*2 


1.48 

0.328 


*10 percent of authority. 


TABLE II. RELATIVE CONTROL EFFECTIVENESS: 
AERODYNAMIC CONTROLS 

[a = 15°; TAS = 133 knots] 

(a) Relative control effectiveness matrix, original controls 


<5o 

Sr 

S D 


0.220 

0.886 

0.409' 

Roll 

0.040 

0.965 

0.258 

Spiral 

0.660 

0.716 

0.229 

Dutch roll 


(b) Control-mixing matrix* 



12.2 

31.4 

-25.1 

14.7 

9.03 

7.57 


*A11 values are given in degrees. 


(c) Relative control effectiveness matrix, pseudo controls 



v 2 


0.994 

0.107' 

Roll 

0.914 

0.407 

Spiral 

0.297 

0.955 

Dutch roll 
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TABLE III. RELATIVE CONTROL EFFECTIVENESS: 
AERODYNAMIC AND THRUST-VECTORING CONTROLS 

[a = 15°; TAS = 133 knots] 


(a) Relative control effectiveness matrix, original controls 



6r 

S D 

5 r,roll 

<5v,yaw 


0.178 

0.719 

0.332 

0.029 

0.583' 

Roll 

0.032 

0.768 

0.205 

0.001 

0.606 

Spiral 

0.583 

0.633 

0.202 

0.074 

0.461 

Dutch roll 


(b) Control-mixing matrix* 


( 6 a 


' 12.1 

31.2- 

Sr 


-20.0 

11.7 

So 

► = 

7.98 

8.07 

^t»,roll 


1.35 

3.07 

^ ^u,yaw • 


.-8.24 

3.84. 


*A11 values are given in degrees. 


(c) Relative control effectiveness matrix, pseudo controls 



V2 


0.988 

0.155' 

Roll 

0.918 

0.398 

Spiral 

0.416 

0.909 

Dutch roll 


TABLE IV. CONTROL DEFLECTIONS: AERODYNAMIC AND THRUST- VECTORING 

[a = 15°; TAS = 133 knots] 


Control 

variable 

(authority) 

Control deflection, degrees 

Baseline 
system 
(fig. 3(a)) 

Aerodynamic 
controls 
(fig. 3(b)) 

Thrust-vectoring 
controls 
(fig. 3(c)) 

Sa (40°) 

*6 

2.19 

2.72 

Sr (30°) 


*-4.5 

*-4.5 

So (20°) 

*3 

1.62 

1.79 

S v , roll (30°) 



0.304 

s v , yaw (15°) 



-1.85 


*15 percent of authority. 
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p, deg/sec 


0 

-12 

100 

4> , deg 

0 

-100 

12 

r, deg/sec 

0 

-12 

3 

3 » deg 
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-3 


— »—| }-* — l second 
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Controls on 


(a) Ailerons and differential tail (lateral controls). 

Figure 1. Response of linear airplane model to baseline controls, a = 10°. 
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p, deg/sec 

0 -- 
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• » deg 


100 T 


0 + 


-100 -*■ 


12 T 

r, deg/sec 

0 ■■ 
-12 -• 


e. deg 

0 -- 
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Controls on 


(b) Rudder (directional control). 
Figure 1. Concluded. 
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1 second 


p, deg/sec 


,, deg 


r, deg/sec 


i, deg 


12 T 


0 + 


-12 x 


100 T 


0 + 


-100 


12 r 


0 + 


-12 •*- 


3 T 


0 + 



-3 J- 

H H 

Controls on 

(a) Pseudo control t>! (lateral control). 

Figure 2. Response of linear airplane model to pseudo controls, a = 10° 
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p, deg/sec 


<|>» deg 


r, deg/sec 


6, deg 


— }-« — 1 second 



Controls on 


Controls on 




(a) Baseline system. 


(b) Aerodynamic controls, (c) Aerodynamic and 

thrust-vectoring 

controls. 


Figure 3. Responses to step lateral-control commands, a = 15°. 
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